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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 

TECHNICAL MEMORANDUM NO. 1231 

REPORT ON INVESTIGATION OF DEVELOPED TURBULENCE* 

By L. Prandtl 


TECH LIBRARY KAFB, NM 



The recent exper ime nts By Jakob and Erk (reference l) on the 
resistance of flowing water In smooth pipes, which are in good agreement 
with earlier measurements By Stanton and Pannell (reference 2) , hare 
caused me to change my opinion that the empirical Blasius law (resist- 
ance proportional to the 7 A power of the mean velocity) was applicable 
up to arbitrarily high Reynolds numbers. According to the new tests 
the exponent approaches 2 with increasing Reynolds number, where it 
remains an open question whether or not a specific finite limiting 
value of the resistance factor X is obtained at R - <». 

With the collapse of Blasius* law the requirements which produced 
the relation that the velocity in the proximity of the wall varied in 
proportion to the 7th root of the wall distance muBt also become void 
(reference 3). However, it is found that the fundamental assumption 
that led to this relationship can be generalized so as to furnish a 
velocity distribution for any empirical resistance law. These funda- 
mental assumptions can be so expressed that for the law of velocity 
distribution in proximity of the wall as well as for that of friction 
at the wall, a form can be found in which the pipe diameter no longer 
occurs, or in other words, that the processes in proximity of a wall 
are not dependent upon the distance of the opposite wall. 

For the velocity u (time average value) at y distance from the 
wall only one nondimens ional can then be formed, namely, (v = kine- 
matic viscosity) , giving for u a formula of the form 

u W 

where C is a velocity and cp an arbitrary function. The shearing 
stress at the wall then must be 

t = £pC 2 (2) 

where £ is a constant. 


*"Bericht iiber Untersuchungen zur ausgebildeten Turbulenz." 
Zeltschrift fur angewandte Mathematik und Mechanlk, vol. no* 2, 
April 1925, pp- 136-139. 
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To obtain the function cp we proceed from the value 
by theory would be constant = 7 and put 


d Zn y 
d Zn u* 


which 


d Zn y 
d Zn u 


= f(o) 


(3) 


where 0 = Zn^L. With Zn u = tj and Zn y = 0 - tj + Znv, formula ( 3 ) 
becomes 


do 

dtj 


f(o) + 1, 


and, after integration and removal of logarithms. 


u 


Ce 


J“ «£_ 

f(o; + 1 


which, after including 



( 4 ) 


C) 


gives the velocity profile in parameter representation. 

The' empirical law for wall friction reads t= with k a 

function of ~ (u = mean velocity, a = pipe diameter) . 


Discounting for simplicity the difference between the mean velocity 
and the velocity in the center up and assuming for it the value from 

formula (4), although it is no longer exactly true in the pipe center, 

u-]_a 

we put with Op = Zn-^— 


~ Ink = g(p ± ) 


By (2) and (4) we get 



r 

2 

pup 




do 

f(a) + 1 


( 6 ) 
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or 


Zn£ 


= “ g(ai) 


+ 2 


' 1 da 

f(a) + 1 


const. 


(T) 


hence after differentiation with respect to cr-j. 


= 1 + f(a x ) ” - 2 g'fa) g) (®) 

With this the problem is solved and it is readily seen that f(cr), 

which for g'(cr) = ^ assumes the value 7 as before, increases with 

decreasing g'(cr). A more accurate experimental check is awaited, but 
even so it is plainly seen that at Reynolds numbers of about 200,000 
the 8th root of the wall distance is definitely better than the 7th root* 
For g'(a) = 1, f(cr) = 1, which corresponds to the laminar boundary 
zone. 


Furthermore I would like to speak of a formula intended for a 
hydrodynamic calculation of the distribution of the base flow of a 
turbulent motion under the most varied conditions. 


After various fruitless attempts gratifying success was attained. 
In addition it was found that the formula for the apparent shearing 
stress r produced by the interchange of momentum, lends itself to a 
very clear explanation. 


In the Boussinesq formula 
turbulent 


t = 


du 

: <3y' 


s is a measure for the 


"exchange" and in its dimension, which is the same as that 
of V, It is the product of a length and a velocity. The velocity is 
the transverse velocity w at which on the average the fluid bodies 
advancing from both sides pass through the layer with the t ime average, 
value of the velocity = u. 


The liquid bodies coming from the side of the greater velocity 
entertain higher values of velocity u, those from the side of smaller 
velocities, smaller values, with the result that more* and more momentum 
Is transported In one direction than in the other (excepting the point 
of Umax) • The desired length Z is characterized by the fact that it 
indicates the distance of the particular layer, in which the average u 
velocities, which the liquid bodies have at their passage, are found as the 
time average value of the flow velocity. Approximated these velocities 
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are u + Z— and u - Z~. (Incidentally Z is in agreement in order 

of magnitude with, the diameter of the fluid bodies (more accurately it 
is the decelerating path of the fluid bodies in the remaining fluid, 
which is, however, proportional to the diameter).) As to the length Z 
it can, for the present, only be stated that it must approach zero at 
the wall, where only bodies of smaller diameter than the w all distance 
can move as discussed. Elsewhere Z is to have a very regular distri- 
bution. If 0 is the average proportional share of the surface 


occupied by the fluid bodies entering from one side, a momentum 


d du 


per second passes at this side through the unit surface, and approxi- 
mately the same amount from the other side. This confirms the Boussinesq. 
theorem, so we can put e = 20wZ. 


The next problem is to find a practical formula for the mi -ring 
speed w. This mixing speed is rapidly reduced and must be 
continuously renewed. Hence the assumption that it is produced at the 
concurrence of two bodies of different velocity u and therefore 
proportional to the velocity difference, that is, the ma gn itude of Z— . 

With this, however, if all unknown numerical factors are thrown on ^ 
the not more accurately known length Z, the apparent shearing stress 
becomes 


T 


pZ 2 


du du 
dy dy 


( 9 ) 


This formula necessitates a correction for the case that — = 0. 

dy 


For the creation of velocity w the neighborhood In a certain 
width cooperates; It does not become zero when ^ = 0} It rather can 


be put proportion al to a . statistical average value of 
v2 

proportional to 


du 

dy 


that Is 


r du \ 

— ) . If the velocity profile varies in flow direc- 


tion, as in convergent divergent channels, the points over which the 
averages are made must be shifted upstream for a certain amount, since 
the process of development of the velocity w takes time* 


Foimiula ( 7 ) has already proved itself in many respects. In a pipe, 
for Instance, the shearing stress is according to the equilibrium con- 
ditions proportional to the distance r from the center, hence 
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Assuming l as constant gives 


u = A - Br 


3/2 


( 10 ) 


for r > 0, the reflection for r < 0 (fig* 1) • At r = 0 the radius 


du 


the previously 

discussed statistical value that can he approximately -written as 


of curvature is then equal to 0} taking instead of — 


teuf + i'.2/&Y 

W / 


the curvature at 'r = 0 is finite* The actual 


velocity distributions (reference 1, p. 21) exhibit the singular behavior 
of the center very plainly (fig* l) * The refinement just mentioned 
needs to be applied only where increased accuracy requirements are 
involved. The total velocity distribution in the pipe is fairly accu- 
rately obtained -when (in the range of the l/T power law) l is put 

proportional to [(a - r) (a + r)] a = pipe radius. 


The foregoing formulas have also been applied to the case of the 
"free turbulence, that is, flows without confining walls, as for 
example, of a fluid jet diffusing in a chamber and to the intermingling 
of a homogeneous air stream with the adjacent air at rest, which 1 b 
entrained by it. For stationary flows of this type the formula l = cx 
is appropriate, x being the distance from the point where the .-mixing 
starts. The case reproduced in figure 2 leads to the differential 
equation for the stream function Ff^ 


2cF"F m ± FF" = 0 

which is solved by F" = 0 or by 2cF MI ± F = 0. The two solutions of 
uniform velocity and variable velocity abut in F 1 " with discontinuity. 
This and other calculations were numerically carried out by Tollmlen, 
who is to publish an article on it- The agreement of his calculations 
with experimental data 1 b excellent. 

Further experimental studies included the velocity distributions 
in channels of other than circular sections with very unusual results 
and for which the explanation has not yet been found* Slightly 


^It is approximately so that u^ plotted against the 
section of the channel gives a sloping surface* 


cross 
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divergent a n d convergent smooth-vailed channels have also been investi- 
gated, rough -vailed channels are in preparation. It is hoped that our 
formula will prove itself here also. 


Translated by J. Vanier 
National Advisory Committee 
for Aeronautics 
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